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1. INTRODUCTION 

Triplet pairing is exemplified physically in at least two strongly interacting quan- 
tum many-fermion systems: 

(i) In laboratories on earth: Superfluid phases are observed when liquid ^He is 
cooled to mK temperatures (triplet pairing of spin-1/2 atoms) [1]. 

(ii) In the quantum fluid interior of a neutron star: It is generally believed that 
the neutrons pair-condense at temperatures below ~ 10^ K (triplet pairing of 
spin-1/2 nucleons) [2]. 

There are important differences in these realizations. In dense neutron matter, 
^Sq pairing, which dominates in the inner crustal region of the star, is quenched by 
the strong short-range repulsion in the nucleon-nucleon [NN) interaction. At A^A^ 
collision energies corresponding to the relevant density regime, the most attractive 
phase shift is found in the channel, indicating that this state will dominate the 
pairing. The ^Pq phase shift is only weakly attractive, while the '^Pi phase shift is 
repulsive; these channels are therefore expected to have only minor effects on the 
problem. The state dependence of the spin-orbit component of the A^A^ interaction 
is responsible for this pre-eminence of the state. A complication is introduced 
by the tensor component of the interaction, which implies a coupling of the ^P2 and 
■^^2 channels. 

On the other hand, in liquid ^He, the atom-atom interaction is spherically sym- 
metric and spin-independent to a very good approximation. Tensor and spin-orbit 
components are effectively absent. Thus there is no coupling to higher L waves, but 



the three J states for L = 1, S' = 1, namely ^P25 ^Pii and ^Pq now enter on an equal 
footing. 

On balance, the neutron matter - in spite of the complication of tensor coupling 
- is simpler than liquid ^He, since fewer magnetic substates enter the problem. 
Therefore we focus first on neutron matter as a "training ground" for the more 
challenging case of liquid ^He. 

Quantitative, microscopic prediction of the superfluid phases and phase diagram 
of systems manifesting triplet pairing (or pairing in orbital angular momentum states 
beyond S'-wave), has hitherto eluded the best efforts of theorists. For one thing, the 
popular Ginzburg-Landau approach is restricted to temperatures near the critical 
temperature Tg. Secondly, iterative procedures, routinely employed to calculate the 
energy gap in weakly interacting systems with S'-wave pairing, are afflicted with 
slow convergence and uncertain accuracy when applied to the many coupled, sin- 
gular nonlinear integral equations that come into play for pairing in higher angular 
momentum states. The limitations of standard iterative procedures become most ap- 
parent when one attempts to construct the superfluid phase diagram of the system, 
which is controlled by tiny energy splittings between different solutions of the BCS 
pairing problem. 

These difficulties can now be overcome by means of a recently introduced sepa- 
ration method [3-8] for solving the BCS gap equations associated with pairing states 
of arbitrary angular- momentum content. In this paper, we implement this method 
to identify and classify the spectrum of phases of the ^P2-^F2 pairing model. Nomi- 
nally, the model refers to superfluid neutron matter at densities of order twice that 
found in heavy nuclei. However, the results obtained are universal in the sense that 
the angular properties of the solutions corresponding to the existing phases are inde- 
pendent of the detailed nature of the system under study, and especially its pairing 
interaction. Indeed, these results can be derived analytically and depend only on 
two energy scales, one governing the overall strength of the pairing effect and the 
other measuring the strength of the coupling between P and F states. There exist 
the usual solutions involving a single value of the magnetic quantum number M (and 
its negative) - obviously three in number corresponding to M = 0, ±1, ±2. In addi- 
tion we find that there exist ten real multicomponent solutions. Five of these have 
angle-dependent order parameters with nodes (and therefore are of relatively high 
energy), while the other five are nodeless. In contrast to the case of superfiuid ^He, 
transitions occur between phases with nodeless order parameters. Results for the 
temperature dependence of the competition between the various phases have been 
obtained. In principle, and presumably also in practice, the same approach can be 
used to complete the catalog of superfiuid phases of liquid ^He. 

2. BCS FORMALISM FOR ANY-CHANNEL PAIRING 

Adopting a spinor representation, the 2x2 gap matrix for the general BCS 
problem of pairing in a state of triplet spin (S = 1) and triplet isospin (T = 1) has 
the expansion 

A«^(p) = J] A^ip) (Gf,(n))^^ (1) 

J,L,M 



in terms of the spin-angle matrices 



(Gf.(n))^^= Y: Cl%C^ii,,M.Y^MM . (2) 

Ms Ml 

The coupled partial- wave gap components A£^(p) solve the set of gap equations 

tanh(S(pi)/2T) , . , x 9, , 
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The quasiparticle energy -E'(p) = [^^(p) -|- D^(p)] is constructed from the gap 
components Aj^-^ [p) through 

D\p)= (Ar(p))*A£^Hp)^if/^^^(n) (4) 

LJMLiJiMi 

together with the single-particle spectrum ^{p) of the normal Fermi liquid, often 
parametrized in terms of an effective mass M*. Angular dependence is introduced 
into the quasiparticle energy E{p) via the spin trace 

SJMJ.M. (n) = TV [(Gi^(n))* G^^^^ (n)] , (5) 

which obviously complicates explicit solution of the system of gap equations. 

The pairing matrix elements are generated by the spin-angle expan- 

sion 

V{p,pi)= (-^)'^^\^LL'\Pi)GZin){Gf!,j{m)y (6) 

LL'JM 

of the totality of vertex diagrams irreducible in the particle-particle channel. 

A close inspection reveals that complete solution of this set of equations presents 
awesome difficulties: note, in particular, the coupling to the generic angular momen- 
tum labels Li, Ji, and Mi quantum numbers via the squared-gap quantity -D^(p) 
appearing in the quasiparticle energy E{p). Practical approximate solution (which 
might still be extremely accurate) will require a series of justifiable simplifications. 

To wit, in dealing with neutron matter, the free-space two-body interaction has 
the salient feature that the components of the central forces nearly cancel each other, 
as reflected in the behavior of the experimental P-scattering phases. It is assumed 
that this feature is preserved by the effective interaction inside neutron matter. We 
shall later comment on the veracity of this assumption. If it holds, the pivotal 
role of the spin-orbit force in promoting the pairing channel then implies that 
contributions to triplet pairing from "nondiagonal" terms with L' , Li 7^ 1 or Ji 7^ 2 
on the right-hand side of the set of gap equations can be evaluated within perturbation 
theory, in terms of the set of principal gap amplitudes Af^ (p), with M — 0, ±1, ±2. 



In fact, if we choose to invoke time-reversal invariance, the problem may be 
treated in terms of only three complex functions, namely Af^ (p) with M — ,1 ,2. 

Another simplifying feature is the existence of a small parameter across the range 
of interesting pairing problems. When dealing with anisotropic gaps as may arise in 
pairing states beyond the S-waye, it is helpful to define the energy gap parameter 
Ap as the square root of the angle average of D^(p), evaluated at the Fermi surface. 
The ratio dp = Ap/ep of the gap to the Fermi energy provides a small parameter for 
BCS pairing theory. For (or ^P2-^F2) pairing in neutron matter, a leading order 
approximation in dp is good to about one part in 10^ or 10^. A similar or higher 
accuracy is to be expected in liquid ^He. 



3. THE SEPARATION METHOD 

The separation approach developed in Refs. [3-8] facilitates essential further 
simplifications (linked in part to the existence of a (very) small parameter). In this 
approach, any given pairing matrix element is expressed identically as a separable 
part, plus a remainder that vanishes when either momentum argument is on the 
Fermi surface. We write 

{p\VIl'\Pi) = viL><i>iL'{p)<t>iL'{Pi) + WiL,{p,p^) (7) 

and verify that the choices -yf^, = ((Pf\VIi^,\pp) and (pLL'iP) = {p\^ll'\Pf) /vLl' 
produce the desired behavior when p or pi hits pp. Substitution of this identity 
into the set of gap equations, followed by simple manipulation and argumentation, 
establishes a decomposition of the form 

Ai^{p)=Di''xiAp) (|M| = 0,1,2) (8) 

for the general gap component and provides separate equations for the shape factor 
XLL'iP) (normalized to unity at p = pp) and the numerical amplitudes Dj^^ that 
determine the angle-dependence of the gap. 

The shape factor Xll'^p) given by a nonsingular integral equation involving 
a set of kernels proportional to the residual interaction W[j^,{p.,pi) (see Ref. [8] for 
details). The vanishing of W at the Fermi surface removes the logarithmic singu- 
larity characteristic of BCS theory, which is banished to the set of equations for the 
amplitudes -D^^, which are coupled nonlinear equations for a set of numbers. 

To very high accuracy, the integral equation for Xll' linear and independent 
of the amplitudes Dj^^ . For all practical purposes (e.g., appealing to extremely rapid 
convergence of small-parameter expansions), we are free to make the replacements 

^(P)^ie(p)|, tanh(£;(p)/2T) ^ 1 (9) 

in any integral involving the residual interaction W as a factor. The shape factors 
XLL'ip) therefore determined by the interaction W and may be calculated in- 
dependently of the D'l^ by a matrix inversion routine. The nonlinear and singular 
aspects of the problem reside entirely in the set of equations for the amplitudes D'l^ , 



which may also be solved numerically by standard methods (e.g. Newton-type algo- 
rithms). In fact, we have even been able to solve scaled versions of the D'l^ equations 
analytically in interesting will emerge below. 

In view of these crucial simplifications revealed and expedited by the separation 
transformation, it is seen that 

(i) For given L, L', and J, the x function is universal, i.e., independent of magnetic 
quantum number M and temperature T. 

(ii) The factorization 

A?^(p) = ^rx(p) (|M| = 0,1,2) (10) 
may be asserted for the principal gap components, where xip) = Xiiip)- 

(iii) Explication of the phase diagram of dense superfluid neutron matter reduces to 
determination of the three amplitudes , since the character of the phase 
diagram itself is independent of the shape factor x{p)- 



4. PERTURBATIVE TREATMENT 

Now let us focus on the nondiagonal contributions to the right-hand side of the 
gap equations (3), where nondiagonal means one or more of the conditions L' = 1, 
Li = 1, Ji = 2 is not met. Two nondiagonal contributions are found to be of leading 
importance: The first contains the integral of the product V^iSl^"^^^ /\\^'^ while 
the second contains the integral of the product ViiSl^"^^^ l^^^ . There is a single 
small factor in each product, namely V^i in the first instance and Ag"^^. All other 
nondiagonal contributions involve two or more small factors. Retaining the leading 
nondiagonal pieces and ignoring all the others, we have what is called the '^P2-^-^2 
pairing problem, which is generally regarded as a satisfactory model of superfiuid 
neutron matter. 

Rapid convergence of the pertinent nondiagonal integrals is instrumental to suc- 
cess of the perturbation approach used to treat the nondiagonal effects. To affirm 
this behavior, we observe that the dominant contributions to these integrals come 
from the vicinity of the Fermi surface. Then, if -E'(p) is significantly larger than the 
energy gap value (as it will be for large p), -E'(p) and |^(p)| will practically coincide, 
and consequently the angular integration yields zero. 

We conclude that in evaluating the nondiagonal contributions it suflBces to know 
the "minor" gap components A'^^ (p) at p = pp- These quantities will in fact be 
determined in terms of the coefficients . To excellent approximation, we may 
retain only the dominant contribution to the right-hand side of the gap equations (3) 
that diverges like hidp as Ap — 0. We are then led to the desired connection 

Ar{p=PF)=vDi^, (11) 

where r] = —{pf\Vi^\pf)/vf and vf = {pf\Vii\pf)- Similar relations may be ob- 
tained for other minor gap components, notably and A^-^^. 

What can we say about the parameter rj, which represents the coupling to the 
"nondiagonal" states? If perturbation theory is to be valid, it should be com- 
fortably small. For pairing matrix corresponding directly to in-vacuum neutron- 
neutron interaction, rj depends smoothly on density, varying around 0.3 in the interval 



Pq < p < 2po- Medium modification of the spin-orbit force is probably not impor- 
tant. Due to the relativistic origin of this component, it should not be much affected 
by polarization or correlation corrections. This claim is consistent with empirical 
analyses of the spin-orbit splitting in finite nuclei. 

On the other hand, medium modification of the tensor force may be more sig- 
nificant, especially as one approaches the density at which pion condensation occurs. 
Thus, the parameter 77 is somewhat uncertain, although it is still expected to be 
rather small. 

To sketch out the superfluid phase diagram of the system, we need to deter- 
mine the key amplitudes (M = 0,1,2) to leading perturbative order in the 
coupling T}. Exploiting the linear connection (11) between the amplitudes and 
the minor-component values A|^(pj?), simple manipulations applied to the coupled 
gap equations (1) aX p = pp yield 

Dl^^v^Y^D^ j j 0(p) '""^^^g^|Pj/'^^ 5,T^^nn)x(p)/cipdn = r^v.VM 

(12) 

with (j){p) = {p\V^i\pf)/vf, Eo{p) = E{p; r] = 0), and 

rM = E Dl""^ I I Hr""' (n) + (n)l '-^^^^^^^p' dpdn . (13) 

We restrict the search for solutions of these equations to those with real coeflBcients 
D^^, reasoning that solutions with complex D amplitudes will lie at energies high 
enough to make them physically irrelevant. 

Inserting the explicit form of 5'^^^^^(n) from Eq. (5) into Eq. (12), we derive 
a system of three equations for the gap value A p and the two ratios 

X^=Df/DfVQ and = Df / dI'' , (14) 

which serve generally to determine the angular dependence of the gap (or alternatively 
its composition with respect to the magnetic quantum number M). In a notation 
and form compatible with the earlier treatment of the pure problem [4,6], these 
basic equations read 

A2 + vf [A2(Jo + -^5) - Ai Ji - J3] = r]VFr2 , (15a) 
Ai + VF [-(A2 + 1) Ji + Ai( Jo -f- 4J5 + 2J3)/4] = 7]VFri , (156) 

1 + VF [-{X2J3 + Ai Ji)/3 + J5] = llVFTo , (15c) 

with 

tanh(Eo(p)/2T) p'^dpdn 



/o=l-3^^ h = 3xz/2, f3 = 3{2x^ + z^-l)/2, h = {l+3z^)/2, (17) 



z — COS 6 , X = sin 6 cos (p , y — sin 9 sin ip . 



(18) 



Upon setting r] = 0, Eqs. (15a)-(15c) collapse to Eqs. (12)-(14) of Ref. [4], which 
refer to the case ki = = and were solved analytically in Refs. [4,6]. 

If we had not used the separation method, we would have been faced at this point 
with a system of coupled singular nonlinear integral equations for a set of functions, 
not numbers. The great advantage gained is transparent. 

5. SEARCHING FOR MULTICOMPONENT SOLUTIONS 

The basic equations (15a)-(15c) possess the familiar one-component solutions 
corresponding to definite |M|, i.e., |M| = 0,1,2. In addition, there exist multicom- 
ponent solutions whose structure and spectrum we seek to establish, via a two-step 
process. We first note that J5, which is the only Jj integral containing a principal 
term going like lii{ep/ Ap), is responsible for the gap magnitude Ap. On the other 
hand, J5 is irrelevant to the phase structure. 

Step 1. Thus, we begin by eliminating J5 from the first pair (15a)-(15b) of basic 
equations and reduce the number of Jj integrals in each of the pair to two: 

(A2 + l)[3Ai(A2 + 1) Jo - 2(A2 - 2Ai + 6) Ji] = 7]Bi , (19a) 
(A2 + 1)[(A? - 4A2) Ji + Ai(A2 + 1) J3] = VB2 , (196) 

Bi = 2Ai(2A2 + 3)r2 - 4(A2 - 3)ri - 6Ai(A2 + 2)ro , 

^2 = -Air2 +4A2ri -3AiA2ro. (20) 

Step 2. Assuming A2 7^ 1, we perform the rotation 

{x, z) — {—t sin 1) + u cos 1), t cos 1!} + u sin '&) (21) 

and choose the angle so as to remove the integral Ji from the "sanitized" first 
pair of basic equations, (19a)-(19b). Specifically, if C = tani? is chosen to obey the 
algebraic equation 

AiC'-(A2-3)C-Ai = 0, (22) 

then substitution of the transformed integrals Ji converts the key pair of equations 
(19a)-(19b) into 

iX2 + l)[AiJo + A2J3] = vBu (23a) 
(A2 + l)[Ai Jo + A2J3] = -277^2 , (236) 

with 

= ^Ai(l + A2)(2 - C^) - ^{Xj - 2Ai + 6)C , 
A2 = -3Ai (1 + A2)C' - (A? - 2Ai + 6)C . (24) 

We observe that the left-hand members of Eqs. (23a) and (23b) are identical! 
It is in fact just this coincidence that leads to the striking universalities of the pure 




-10 -5 5 10 



Figure 1. Solutions for the coefficient ratios Ai and A2 in the case rj = 0, corre- 
sponding to the uncoupled pairing problem solved in Refs. [4,6]. The particular 
(point) solution (Ai = 0, A2 = 3) is indicated by the solid dot and the degenerate 
solutions by the solid curves. 

pairing problem first revealed in Ref. [4]. Of course, in this problem, which 
corresponds to = 0, we have the special circumstance that the right-hand side of 
each equation vanishes identically, so that Eqs. (23a) and (23b) coincide and yield 
only the single constraint 

(A2 + l)[^lJ0+^2J3] =0 (25) 

on the parameters Ai and A2. 

It ensues that the solutions of the pure pairing problem display remarkable 
universalities as expressed in two kinds of degeneracies, independently of temperature, 
density, and details of the in-medium interaction: 

• Energetic degeneracy. There exists an upper group of states, degenerate in en- 
ergy, whose angle-dependent order parameters have nodes, and a lower group 
with nodeless order parameters. Relative to the absolute pairing energy, the 
splitting between upper and lower states is small (of order 2% in the neutron- 
matter problem). 

• Parametric degeneracy. The multicomponent solutions, which satisfy the spec- 
tral condition [4,6] 

(A^ + 2 - 2A2) (A? - 2A^ - 6A2) = , (26) 
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Figure 2. Multicomponent solutions of the pairing problem, defined by the 
coefficient ratios Ai and A2. Solutions with nodeless (respectively, node-bearing) 
order parameters in the case of pure pairing are indicated by open (respectively, 
filled) circles. 

display a parametric degeneracy with respect to the parameters Ai and A2. Thus, 
as seen in Fig. 1 (which was constructed analytically), the multicomponent so- 
lutions generally define curves rather than points in the (Ai, A2) plane. 

6. BREAKING THE PARAMETRIC DEGENERACY 

The strong parametric degeneracy intrinsic to the uncoupled "^^2 pairing problem 
is lifted in the case of ^P2-^F2 pairing, i.e., at 17 7^ 0. With the coupling turned on, 
the true solutions of the problem are represented by a set of isolated points in the 
(Ai, A2) plane. 

From the analytical standpoint, the underlying "mechanism" runs as follows. 
Upon equating the right-hand members of the double-sanitized pair of basic equations 
(23a)-(23b) in the small-|?7| limit (|?7| infinitesimal), one obtains an additional relation 
between the parameters Xi{r] = 0) and A2(?7 = 0): 



Air2 - (A2 - 3)ri - SAiTq = . 



(27) 



This relation supplements the spectral condition (26) nontrivially and removes the 
parametric degeneracy. The system formed by (27) and (26) is solved by applying 
the same rotation in x — z coordinates as introduced previously. After some algebra, 
one may then obtain the full set of solutions of the coupled-channel ^P2-^-p2 pairing 
problem. A salient feature of these solutions, made explicit via the separation scheme, 
is their virtually complete independence of the temperature T. 

We summarize the results of this analysis with a catalog of the possible solutions 
for "^P2^"^-^2 pairing described in the BCS framework. One simplification, already 
evident in Fig. 1 but also quite general, is that the relevant pairing energies are 
independent of the sign of Ai, so that we only need to consider Ai > 0. The modified 
picture in the (Ai, A2) plane is displayed in Fig. 2. 

First of all, there remain the three well-known single-component solutions with 
\M\ = 0, 1, or 2. Beyond these, the collection of unitary solutions of the ^^2^ 

pairing problem contains ten multicomponent solutions, corresponding to more 
complicated superfluid phases. 

Five of these additional solutions, denoted Ok (k = 1,...,5), have nodeless order 
parameters and include: 

(a) Two two-component solutions 0±3, identical to those found in the pure 
pairing problem, with Ai = and A2 = ±3. 

(b) Three three-component solutions: 

o Two of them, Oi and O4, are associated with the upper branch of A^ — 
2A2 — 6A2 = and have A2 = 3(\/2T— 4)/5 and A2 = 3, respectively. 

o The third, O2, is associated with the lower branch of the same equation and 
has A2 = -3(v^ + 4)/5. 

The other five solutions, denoted by Xk (A; = 1, . . . , 5) do possess nodes. 

(a) Two two-component solutions X±i, again identical to those found in the pure 

pairing problem, with Ai = and A2 = ±1. 

(b) Three three-component solutions, X2, X^, and X4, associated with the parabola 
A2 = Af/2 + 1 and having A2 = 13 - 2^35, A2 = 3, and A2 = 13 + 2^35, 
respectively. 

These general features of the spectrum of solutions of the ^P2-^F2 problem are 
expected to persist even if |?7| is not so small. 



To complete the phase diagram of superfluid neutron matter, we need the gap 
values Ai? for the various phases. These are found through the third basic equation 
(15c), which involves J5. To determine which phase wins the competition at a given 
temperature T, we compare the free-energy shifts 



7. THE BATTLE BETWEEN PHASES 




(25) 



due to pairing in the corresponding superfluid states, where g is the relevant coupling 
constant. 
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Figure 3. Temperature dependence of the splitting between the phase Oi (or 
O2) and the phases 0±z (or the one-component phase with M = 0), in terms of 
the difference SA'^{T) between the respective A|,(T) values, measured relative to 
|5A|(T = 0)|. To set the energy scales involved, the latter quantity is around 2% 
of A'jp for rj — 0.3, while the gap parameter Af{T — 0) itself reaches a maximum 
value near 0.44 MeV at kp — 2.1 fm"-*^ when the pairing interaction is given by the 
Argonne vig potential and free normal-state single-particle energies are employed [8] . 



At low T, the only viable contestants arc the solutions with nodeless order 
parameters, since the other solutions lie too high in energy. The separation between 
the two groups of states simply cannot be bridged if the value of |?7| remains rather 
small. 

As documented in the last section, the parametric degeneracy of the pairing 
problem in the Ai — A2 plane, embodied in the spectral relation (26), is completely 
eradicated when the 77-coupling is switched on. On the other hand, the energetic 
degeneracy between the different superfluid phases is only partially lifted. Specifically, 
the spectrum of pairing energies decays into several groups of nearly degenerate 
states: The Oi and O2 phases form the lowest-energy group, followed by the phases 
0±s along with the one-component phase with M = 0, and so on. 

Raising the temperature from T = toward the critical temperature Tc, the 
splitting between the two groups lowest in energy shrinks, until, at T ~ 0.7 T^, their 
roles are interchanged and transitions occur. This behavior is quantified in Fig. 3. 
We observe that the analogy with the A-B phase transition in superfiuid ^He is 
imperfect. Whereas the order parameter of the A- phase solution has nodes while 
that of the B-phase does not, the transition in neutron matter takes place between 
phases with nodeless order parameters. 




8. CONCLUSIONS 



We close with some details and caveats and examine some prospects for observ- 
able implications of our findings. 

First, there is the matter of unresolvable fine structure in the phase diagram. 
Introduction of the ^Pq and pairing channels would entail a further lifting of 
energy degeneracies, but current ignorance of the in-medium interaction precludes 
reliable estimation of these effects. 

Secondly, we point out that the region near the critical temperature Tc may 
admit subtle phase transitions not uncovered by our analysis. In this regime, strong- 
coupling corrections [9,10] can no longer be neglected in deciding the contest between 
phases. Also, the character of the phase diagram may be influenced significantly by 
external magnetic fields, a fact not considered in our treatment. And further, the 
phenomenon of fermion condensation [11], occurring as a precursor to pion conden- 
sation [12-14], may provide an exotic source of phase transitions in superfluid neuron 
matter. 

In summary, it has been established that the superfluid phase diagram of dense 
neutron matter can in principle exhibit several triplet superfluid phases (at least 13!). 
Transitions between the different phases are expected to occur as a young neutron 
star cools. Since the gap value changes in these transitions, their occurrence may 
ultimately be detected in the thermal history and/or the rotational dynamics of the 
star. Such transitions may produce a signiflcant variation of the moment of inertia 
of the star, or alterations of the distribution of the angular momentum between the 
crust and the vortex system, resulting in a change of the star's angular velocity. 

The existing analysis of triplet pairing in neutron matter provides a firm foothold 
for the more ambitious project of a complete characterization of the repertoire of su- 
perfluid phases of liquid ^He, which is even more complex, promising a richness of 
phenomena not yet revealed by experiment or theory. Since ^P2, ^Pi, and ^Pq states 
all enter the picture, one will be faced with nine nonlinear equations for the relevant 
gap components (in contrast to the five encountered in the pure ^P2 case). Never- 
theless, the separation method will continue to provide an incisive tool of analysis. 

Thus, the procedure applied here may in principle be extended to count all possi- 
ble solutions of the set of nine BCS equations for the gap function of superfluid "^He. 
We may begin by severing the couplings between ^P2, ^Pi, and ^Pq components, 
which facilitates a reference analytic solution. With this solution at hand, we may 
restore the channel couplings and solve the resulting set of equations numerically step 
by step, under increase of the coupling parameters. Attention then turns to (i) veri- 
flcation of the structure of the A-phase predicted by ABM [15,16], (ii) evaluation of 
the difference between free energies of A and B phases vs. temperature and pressure, 
and the barrier between the two phases (important for understanding the huge delay 
of the phase transition between the supercooled A-phase and the B-phasc [17-19]), 
(iii) investigation of strong- coupling corrections beyond the Ginzburg-Landau ap- 
proximation, and (iv) examination of the phase diagram of liquid "^He in an external 
magnetic field. 
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